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Abstract 

In Extensive Air Shower (EAS) experiments, the primary direction is reconstructed 
by the space-time pattern of secondary particles. Thus the equalization of the transit 
time of signals coming from different parts of the detector is crucial in order to get 
the best angular resolution and pointing accuracy allowed by the detector. In this 
paper an off-line calibration method is proposed and studied by means of proper 
simulations. It allows to calibrate the array repeatedly just using the collected data 
without disturbing the standard acquisition. The calibration method is based on 
the definition of a Characteristic Plane introduced to analyze the effects of the 
time systematic offsets, such as the quasi-sinusoidal modulation on azimuth angle 
distribution. This calibration procedure works also when a pre-modulation on the 
primary azimuthal distribution is present. 
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1 Introduction 



In EAS experiments, the space-time information of the secondary particles is 
used to reconstruct the primary direction [1,2,3]. The space information refers 
to the detector unit position while the time information is achieved usually 
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by TDC (Time to Digital Converter). The former is easy to measure and sta- 
ble in a long period, while the latter depends on detector conditions, cables, 
electronics, etc, and usually varies with time and environment. The time off- 
sets are the systematic time differences between detector units, which lead to 
worse angular resolution, and more seriously, to wrong reconstruction of the 
primary direction. As a consequence the azimuthal distribution is deformed 
according to a quasi-sinusoidal modulation [4]. Thus the correction of these 
systematic time offsets [5] is crucial for the primary direction reconstruction, 
much more when the EAS detector is devoted to gamma ray astronomy and 
the pointing accuracy is required in order to associate the signals with astro- 
physical sources. Usually manual absolute calibration by means of a moving 
probe detector is used in EAS arrays, but this method takes time and man- 
power. The difficulty increases taking into account that periodical checks are 
necessary to correct possible time-drift of the detector units due to change in 
the operation conditions. Furthermore the number of detector units in current 
EAS arrays is getting larger and larger. As a conclusion, effective off-line cali- 
bration procedures are greatly needed because they do not hamper the normal 
data taking and can be easily repeated to monitor the detector stability. 

Here a new off-line calibration procedure is presented. It does not depend on 
simulation and is very simple in the case of a uniform azimuthal distribution. 
It works also when some small modulation of the azimuthal distribution is 
expected, for istance due to the geomagnetic field. The correctness of this 
calibration method has been checked by means of simple simulations both in 
the case of uniform and modulated azimuthal distribution. 



2 Characteristic Plane 

In EAS experiments, for an event % the time Uj is measured on each fired 
detector unit j, whose position (xj, i/j) is well known. The primary direction 
cosines /j = sin^cos^, rrii = sin#jsin0j (0, and fa are zenith and azimuth 
angles) can be reconstructed by a least squares fit. Taking into account the 
time offset Aj typical of the detector unit and assuming that the shower front is 
plane and the time-spread due to its thickness is negligible, the plane-equation 
is 

c(tij - Aj -toi) = kxj + miUj (1) 

where c is the light velocity, and toi is another parameter of the fit. But the 
time offset Aj is unknown and the goal of the calibration is just to determine 
it. A traditional off-line calibration method is based on the study of the time- 
residuals but their removal does not guarantee the removal of the complete 
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offset. Therefore one can assume that the time offset Aj is the sum of two 
terms: the residual term and another unknown term. Being unaware of Aj the 
plane-equation goes like: 

c (Uj - 4) = l 'i x 3 + m'iVj ( 2 ) 



giving the fake direction cosines l\ = sin 9[ cos 0^, vn! i = sin 0[ sin 0-. From Eq.s 1 
and 2 and neglecting the residuals, it results: 

Aj = op- + bp + 5 0l (3) 

c c 



where aj = l[ — k, bi — w! i — rrii, and 5oi = t' Qi — t 0i is an irrelevant time- shift 
equal for all the units. One can conclude that the offset Aj is correlated with 
the position of the detector unit. The quantities a i; bi define a Characteristic 
Plane (CP) in the (x, y, A) space, depending only on the fired unit pattern, 
representing the difference between the reconstructed plane without consid- 
ering the time offset (FP: Fake Plane) and the real one (RP: Real Plane). 
Events firing different sets of units have different CPs, while events firing the 
same set of units have the same CP, that is the difference between the FP and 
RP is the same. We define the CP of an EAS array like the average difference 
between FPs and RPs, i.e. the systematic deviation between FP and RP (the 
pointing accuracy). The CP is fully determined by the direction cosines 

a = (I') — (I) = sin cos 0o, b = (mf) — (m) = sin# o sin0 o (4) 



associated to the angles 9 , O - 



2. 1 Quasi- Sinusoidal Modulation 



If the probability density function (PDF) of the primary azimuth angle is 
f((f>\0), one can deduce that the presence of the CP introduces a quasi-sinusoidal 
modulation of the reconstructed azimuth angle distribution : 



fW) = me) 



1 + 



r 2 sm 2 ((b' 



cos(0' - 00 ) 



(5) 



0oJ 



where r = sin^ /si n ^- The PDF of the reconstructed azimuth angle is a 
combination of multi-harmonics of odd orders with the the amplitude approx- 
imately proportional to r 2n+1 (n = 0,1,2...) when r << I. The time offset 
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does not introduce even order modulations into the reconstructed azimuth an- 
gle distribution. When f{4>\0) = 1/2-7T, the first harmonic becomes dominant 
and the PDF of the reconstructed azimuth angle goes as 

f(^) = ^[l + rcos(0'-0 o )] (6) 



One can observe that the modulation parameters depend on the angles 9q and 
(f> connected to the CP (see Eq.s 4). The phase is just 4> , while the amplitude 
is proportional to r. By integrating f'(<f)'\6) over 9 it results 



1 

2^ 



1 + sin 9 



1 

sin 9 



cos 



(7) 



A fast Monte Carlo simulation was done to check the above conclusion. The 
azimuth angle was sampled uniformly over [0, 27r] and the zenith angle from 
a typical distribution modulated according to cos 6 9 (the mode value is ~ 22° 
and (1/ sin#) = 3.44). CPs with different 9 and 0o were assumed, subtracting 
sin 9q cos 0o and sin 9q sin 0o from the original direction cosines, respectively, in 
order to get the new direction cosines. Fig. 1 shows the reconstructed azimuth 
distributions for two different CPs with 9q = 0.02 rad, <fto = 0.4 rad and 9 = 
0.20 rad, <fio = 1.2 rad, respectively. The first distribution is well reproduced by 
a best-fit function like that of the Eq. 7 as expected for small values of 9q. Also 
the fit parameters are in agreement with the simulation parameters. The higher 
order harmonics must be taken into account in order to well reproduce the 
second distribution (A2 and A3 are the amplitudes of 2nd and 3rd harmonics) 
because in this case 9q and r are larger. 



in 

c 180L 



xirr 



CP e : 0.02, %: 0.4 


CP 6 :0.20, <t> ; 1.2 


Entries 1e+07 
X 2 /ndf 117.4/97 
Constant 1e+05±32 
% 0.01996 ±0.00013 
*„ 0.4009 ± 0.0065 


Entries 1e+07 
X 2 /ndf 102.5/95 
Constant 1e+05±32 
8 0.1 938 ±0.0001 
1.201 ±0.001 
A2 0.1 336 ±0.0004 
A3 0.01 899 ±0.00042 






Fig. 1. Azimuth angle distributions fitted with harmonic functions (see the text for 
comments). 
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3 Characteristic Plane Method 



According to Eq. 1, if Zj and m 8 were exactly known, then any event can be 
used to relatively calibrate all the detector units hit by that shower, while 
l[ — a and m^ — b can be taken as unbiased estimate of U and m^. Therefore 
the time correction is determined by a and b, i.e. the CP of the EAS array, 
according to Eq. 3. 

Suppose that the primary azimuth angle is independent on the zenith angle 
and distributes uniformly, then (I) = 0, (m) = 0. Thus a = (I') — (I) = (l f ), 
b = (ml) — (m) = (ml), which means that the CP of an EAS array can be 
determined by the mean values of the reconstructed direction cosines. Then 
the time offsets can be calculated by means of the off-line analysis of the 
collected data. 

3.1 A simple simulation as a check of the CP method 

Another fast geometrical simulation was implemented in order to check the CP 
method. One million of showers were extracted from the same distributions 
of 6 and (f> used in Sec. 2.1. The arrival primary directions were reconstructed 
by an array of detector units (10 x 10 units on a surface of 40 x 40 m 2 ). The 
times measured by each unit were shifted by systematic time offsets (first plot 
of Fig. 2). As a consequence the primary directions were reconstructed with 
respect to a CP with a = 3.97 x 10~ 3 and b = 7.34 x 10~ 3 (mean values of 
the reconstructed direction cosines). From the Eq.s 4 it is trivial to estimate 
9 = 8.3 x 10~ 3 rad and O = 1-07 rad. 

The reconstructed azimuth distribution is fitted according to Eq. 7 (see the 
first plot of Fig. 3). As expected the modulation coefficient pi and the phase p2 
are compatible with sin#o(l/ sin#) = 0.029 and 0o> respectively. The angles 
between reconstructed and "true" direction are shown in the first plot of Fig. 4. 

The calibration based on the CP method allows to correct the time measure- 
ments, removing the effect of the time offset on each detector unit. In the 
second plot of Fig. 2 the offset- calibration differences are almost null and the 
RMS is lower than 2.3 x 10~ 2 ns. As an effect of the CP calibration the mod- 
ulation disappears in the azimuth distribution (second plot of Fig. 3) and the 
reconstructed directions are very close to the "true" ones (see the second plot 
of Fig. 4). Then the validity of the CP method is fully confirmed. 
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Fig. 2. First plot: systematic time offsets introduced in the simulation of the time 
measurement. Second plot: differences between systematic offset and calibration 
correction. 
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Fig. 3. Azimuth distribution (first plot: before CP calibration, second plot: after CP 
calibration). In the first plot the sinusoidal fit is superimposed. 




Fig. 4. Angles between "true" and reconstructed directions (first plot: before CP 
calibration, second plot: after CP calibration). 
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3.2 P re-modulation on the primary azimuth angle 



The assumption for the CP method is that the mean values of the primary 
direction cosines are null. Generally this is not true for EAS experiments. The 
possible primary anisotropy, the detection efficiency depending on the azimuth 
angle, the geomagnetic effect, and so on, introduce pre-modulation into the 
azimuth angle distribution. Assuming that the 0-distribution is independent 
on 9, the pre-modulation can be described typically as: 



1 + H 9n COs(n0 + 0„) 
n=l 



(8) 



Only gicos(0 + 0i) contributes to the mean values of the primary direction 
cosines. Therefore they result 

(I) = cos 0i (sin#) , (m) = — ^ sin0i (sin0) . (9) 



The CP method annulls (I) and (m) leaving a sinusoidal modulation on the 
distribution of the new 0" azimuth angle. When g\ and gi are small enough 
and the higher order harmonics can be ignored (see Sec. 2.1) the distribution 
approximately is 

f (0") = 7^ [1 + 9'i cos (0" + 0l) + g 2 cos (20" + 2 )] (10) 



where 



9i = 9i 



2 X 'Xsin^ 



On the basis of this result one can conclude that the calibration with the 
CP method does not remove completely the pre-modulation on the primary 
azimuthal distribution. The g±, #2 amplitudes and the 0i, 02 phases can be 
determined from the reconstructed azimuth angle distribution according to 
Eq.s 10 and 11. Then the direction cosines of the real CP can be determined 
by subtracting the pre-modulation term (Eq.s 9). 

Fast simulations have been used also to check the calibration method in the 
case of pre-modulation with one and two harmonics (gi = 0.05, 0i = 0.3 rad 
and g 2 = 0.02, 2 = 1.2 rad). The results are very similar to those of Sec. 3.1 
confirming that the method works also when a pre-modulation is present. In 
Fig. 5 the "true" azimuthal distribution and the distribution after the first 
step of the calibration are shown. As expected the second distribution is well 
reproduced by Eq.s 10 and 11. 



7 




Fig. 5. The pre- modulation "true" azimuth angle distribution is fitted by the 
function /(</>) = k [1 + g\ cos(0 + 4>\) + 52 cos(20 + 4>2)\- The after-calibration az- 
imuthal distribution is fitted according to Eq.s 10 and 11. The fit-parameter values 
are in full agreement. 

3.3 Geomagnetic Effect 



The geomagnetic field inflects the charged primaries and leads to the well 
known East- West effect (with the modulation period of tt which does not 
modify the mean values of the reconstructed direction cosines and does not 
invalidate the CP method), while the secondary charged particles of EAS 
are separated in the geomagnetic field with the lateral distribution getting 
wider and flatter, thus affecting the detection efficiency [6]. A non- vertical 
geomagnetic field destroys the uniformity of the detection efficiency along the 
azimuth angle which will further leads to quasi-sinusoidal modulation on the 
azimuth angle distribution [7]. The geomagnetic effect on the secondaries is 
typically the most significant pre-modulation (as described in Sec. 3.2) with 
amplitude of the order of few percent and very slight variations with the zenith 
angle. This is just the case discussed in the above section and the modulation 
can be determined according to Eq. 10, after which the time can be off-line 
calibrated using the CP method. 



4 Conclusion 



The definition of the CP makes it easier to understand the effects of the de- 
tector time offsets in EAS experiments, and makes the off-line calibration 
possible. One can successfully correct the time offsets and remove the quasi- 
sinusoidal azimuthal modulation (with <fr depending on 6). The calibration 
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procedure has been analytically defined and checked by means of fast simu- 
lations. The CP calibration is very simple when the "true" azimuthal distri- 
bution is uniform (this feature can be also achieved by selecting special event 
sample). The CP method works also when a "true" pre- modulation (with (f> 
independent on 6) of the azimuth distribution is present. The improvement in 
the pointing accuracy is well shown in Fig. 4 for the simulation of Sec. 3.1. 
In real cases, the pointing accuracy will depend on the detector performances 
and on quality and statistics of the data used for the CP calibration. 

This method has been successfully applied to calibrate EAS detectors [8]. It 
has been also checked [9] by Monte Carlo full simulation and by sampling 
manual calibration. The experimental results of the CP method application 
will be the topic of a future paper. 
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